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Abstract 

We (1) determine the number of Latin rectangles with 11 columns and each 
possible number of rows, including the Latin squares of order 11, (2) answer some 
questions of Alter by showing that the number of reduced Latin squares of order n 
is divisible by /! where / is a particular integer close to \n, (3) provide a formula 
for the number of Latin squares in terms of permanents of (+1, — l)-matrices, (4) 
find the extremal values for the number of 1-factorisations of £;-regular bipartite 
graphs on 2n vertices whenever 1^/c^n^ll, (5) show that the proportion of 
Latin squares with a non-trivial symmetry group tends quickly to zero as the order 
increases. 

1 Introduction 

For l^k^n, akxn Latin rectangle is a k x n array L = (£ij) with entries from 
{1,2, ... ,n} such that the entries in each row and in each column are distinct. Of course, 
L is a Latin square if k = n. We say that L is reduced if the first row is (1,2,..., n) and 
the first column is (1,2,..., k) T . If Rk, n is the number of reduced k x n Latin rectangles 
then Lk,n, the total number of k x n Latin rectangles, is n\ (n— 1)! Rk in /{n—k)\. We will 
sometimes write R ntTl as R n and L n>n as L n . 

The determination of Rk, n , especially in the case k = n, has been a popular pursuit 
for a long time. The number of reduced squares up to order 5 was known to Euler [5] 
and Cayley [I]. McMahon [8] used a different method to find the same numbers, but 
obtained the wrong value for order 5. The number of reduced squares of order 6 was 
found by Frolov [6] and later by Tarry [18]. Frolov [6 J also gave an incorrect count of 
reduced squares of order 7. Norton [H] enumerated the Latin squares of order 7 but 
incompletely; this was completed by Sade [15] and Saxena ^6\. The number of reduced 
squares of order 8 was found by Wells [20j, of order 9 by Bammel and Rothstein [2]. 

The value of Rio was found first in 1990 by the amateur mathematician Eric Rogoyski 
working on his home computer and in the following year by the present first author. The 
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resulting joint paper [12] also presented the number of Latin rectangles with up to 10 
columns. Before he died in 2002, Rogoyski worked for several years on the squares of 
order 11 but the computing power available to him was inadequate, despite his approach 
being sound. Given the advance in computers since then, we can now complete the 
computations moderately easily. 

Several explicit formulas for general n are in the literature (p2], f° r example). Sax- 
ena [TB] succeeded in using such a formula to compute Rj. We will give another very 
simple formula in Section [51 At the time of writing, not even the asymptotic value of 
R n is known. In the case of rectangles, the best asymptotic result is for k = o(n 6 / 7 ), by 
Godsil and McKay [7j. 

2 Terminology 

It can be useful to think of a Latin square of order n as a set of n 2 triples of the form 
(row, column, symbol). For each Latin square there are six conjugate squares obtained 
by uniformly permuting the coordinates in each of its triples. For example, the transpose 
of L is obtained by swapping the row and column coordinates in each triple. 

An isotopism of a Latin square L is a permutation of its rows, permutation of its 
columns and permutation of its symbols. The resulting square is said to be isotopic to 
L and the set of all squares isotopic to L is called an isotopy class. In the special case 
when the same permutation is applied to the rows, columns and symbols we say that the 
isotopism is an isomorphism. An isotopism that maps L to itself is called an autotopism 
of L and any autotopism that is an isomorphism is called an automorphism. The main 
class of L is the set of squares which are isotopic to some conjugate of L. Latin squares 
belonging to the same main class are said to be paratopic and a map which combines 
an isotopism with conjugation is called a paratopism. A paratopism which maps a Latin 
square to itself is called an autoparatopism of the square. 

The number of isomorphism classes, isotopy classes and main classes has been de- 
termined by McKay, Meynert and Myrvold [TTJ for n ^ 10. Our computation does not 
allow us to extract this information for n = 11. However, we do show in Section [7] that 
L n j (6n! 3 ) provides an increasingly accurate estimate of the number of main classes as n 
grows. 

3 The Algorithm 

Our approach is essentially that introduced by Sade [15], adapted to the computer by 
Wells [20l [21], and slightly improved by Bammel and Rothstein [2]. It was also used by 
McKay and Rogoyski [12]. Given a k x n Latin rectangle L, we can define a bipartite 
graph B(L) with vertices C U S, where C = {ci, C2, . . . , c n } represents the columns of L 
and S = {si, s 2 , ■ ■ ■ , s n } represents the symbols. There is an edge from c, to Sj if and 
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only if the symbol j appears in column i of L. Thus B(L) is regular of degree k. Clearly 
B(L) does not determine L in general, since it does not record the order of the symbols 
in each column. For us this is an advantage, since it means there are many fewer graphs 
than there are Latin rectangles. 

Given a regular bipartite graph B on C U S of degree k, let m(B) be its number of 
1-factorizations, counted without regard to the order of the factors. Obviously m(B) is an 
invariant of the isomorphism class of B. In speaking of isomorphisms and automorphisms 
of such bipartite graphs, we will admit the possibilities that C and S are preserved setwise 
or that they are exchanged. (More complex mixings of C and S would, in principle, be 
possible in the case of disconnected graphs, but we have chosen to disallow them.) Using 
this convention, let Aut(S) be the automorphism group of B and let B(k,n) be a set 
consisting of one representative of the isomorphism classes of bipartite graphs B on C U S 
of degree k. 

The theoretical basis of our approach is summarized in the following theorem. Parts 1 
and 3 were proved in [12] and part 2 can be proved along similar lines. 

Theorem 1 

1. The number of reduced k x n Latin rectangles is given by 

R Kn = 2nk\(n-k)\ ^ m(S)|Aut(S)| _1 . 

BeB(k,n) 

2. The number of reduced Latin squares of order n is given by 

R n = 2nk\{n-k)\ ^ w(5)m(-B)|Aut(5)|~\ 

BeB{k,n) 

where B is the bipartite complement (the complement in K n , n ) of B and k is any integer 
in the range ^ k ^ n. 

3. Let B G B(k,n) for k ^ 1. Let e be an arbitrary edge of B. Then 

rn (B) = J2m(B-F), 

F 

where the sum is over all 1-f actors F of B that include e. 

For each k — 1, 2, . . . , 11 in turn we found m(B) for all B(k, 11) using Theorem [U^3) 
and were then able to deduce Rk,n from Theorem [U^l). The number of graphs in B(k, 11) 
is 1, 14, 4196, 2806508 and 78322916, for k = 1, . . . , 5, respectively. For k ^ 6 the 
graphs in B(k, 11) are the bipartite complements of those in — k, 11). The main 
practical difficulty was the efficient management of the fairly large amount of data. Two 
implementations were written in a way that made them independent in all substantial 
aspects (except for their reliance on nauty [10J to recognise the isomorphism class of some 
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graphs). For example, they used different edges e in applying Theorem [D(3), so that gen- 
erally different subgraphs were encountered. The execution time of each implementation 
was about 2 years (corrected to 1 GHz Pentium III), but they would have completed in 
under 2 months if about 3 GB memory had been available on the machines used. 

We also ran the computations for n ^ 10 and obtained the same results as reported 
in [12]. We repeat those results, and include the new results, in Table [TJ It is unlikely that 
R12 will be computable by the same method for some time, since the number of regular 
bipartite graphs of order 24 and degree 6 is more than 10 11 . 

Note that our value of Rn agrees precisely with the numerical estimate given in |12j . 
where estimates of R n were given for 11 ^ n ^ 15. 

4 Some divisibility properties of R n 

Despite obtaining the same value repeatedly for Rn by applying Theorem[T](2) for different 
k in two independent computations, we sought to check our answer further by determining 
its value modulo some small prime powers. By means of the algorithms described in |llj . 
we computed representatives L of all the isotopy classes of Latin squares of order 11 for 
which the order of the autotopism group Is(L) is divisible by 5, 7, or 11. The numbers 
of such isotopy classes are listed in Table [2j Since the number of reduced squares in the 
isotopy class of L is nn\/\ls(L)\, these counts imply that Rn equals 8515 modulo 21175, 
in agreement with our computations. 

We also have the following simple divisibility properties. 

Theorem 2 For each integer n ^ 1, 

1. i?2n+i is divisible by gcd(n! (n— 1)1 Rn, (n+1)!). 

2. R 2n is divisible by n! . 

Proof. Consider R-m+i first. We define an equivalence relation on reduced Latin squares 
of order 2n + 1 such that each equivalence class has size either n! (n— 1)\ R n or (n+1)! . 
Let A be the leading principal minor of L = {l^) of order n. 

If A is a (reduced) Latin subsquare, then the squares equivalent to L are those ob- 
tainable by possibly replacing A by another reduced subsquare, permuting the n par- 
tial rows (£i,n+i, A>+2> • • • >A,2n+i) for 1 ^ K ™, permuting the n — 1 partial columns 
(£ n+ ij,£ n+ 2j, • • • , C-2n+i,j) for 2 ^ j ^ n then permuting columns n + 1, n + 2, . . . , 2n + 1 
to put the first row into natural order. These n! (n— 1)1 Rn operations are closed under 
composition and give different reduced Latin squares, so each equivalence class has size 
n\ (n-l)\R n . 

If A is not a Latin subsquare, the squares equivalent to L are those obtainable by 
applying one of the (n+1)! isomorphisms in which the underlying permutation fixes each 
of the points 1,2, ... ,n. No isomorphism of this form can be an automorphism of a 
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Table 1: Reduced Latin rectangles 
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isotopy classes 
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Table 2: Isotopy classes with certain group sizes 



square in which A is not a subsquare (see [11] Theorem 1]). Hence the squares obtained 
are different and the equivalence class has (n+1)! elements. 

The case of Rm is the same except the second argument gives n\ instead of (n+1)! . □ 

Corollary 1 If n = 2p — 1 for some prime p, then R n is divisible by \_(n — 1)/2J! . 
Otherwise, R n is divisible by |_(n + 1)/2J ! . 

Proof. This follows from Tabled] for n ^ 8. For n ^ 9, note that m \ (m— 2)! for m > 4 
unless m is prime. □ 

Note that, for n ^ 12, the corollary gives the best divisor that can be inferred from 
Table [Hand Theorem [2], except that R13 is divisible by 7! and not merely by 6!. 

Alter p] (see also Mullen [9]) asked whether an increasing power of two divides R n as 
n increases and whether R n is divisible by 3 for all n ^ 6. Theorem [2] answers both these 
questions in the affirmative. Indeed it shows much more — that for any integer m > 1 
the power of m dividing R n grows at least linearly in n. That is, for each m there exists 
A = A(m) > such that R n is divisible by mL A ™l for all n. 

Alter also asked for the highest power of two dividing R n , and here we must admit 
our ignorance. It seems from the evidence in Table [3] that the power grows faster than 
linearly, but we were unable to prove this. 



The literature contains quite a few exact formulas for R n , but none of them appear very 
efficient for explicit computation (though Saxena [IE] managed to compute R-? using such 
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Prime factorisation of R. 



2 
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10 
11 



2 35 • 3 4 • 5 • 2801 • 2206499 • 62368028479 



2 6 • 3 • 7 2 
2 10 ■ 3 ■ 5 • 1103 
2 17 • 3 • 1361291 
2 21 • 3 2 • 5231 • 3824477 
2 28 • 3 2 • 5 • 31 • 37 • 547135293937 



1 
1 

2 2 
2 3 -7 



Table 3: Prime factorisations of R„ for n ^ 11. 



Perhaps the simplest formulas are those in [17] , which relate R n to the permanents of 
all 0-1 matrices of order n. Here we give one that is very similar but uses ±1 matrices 
instead. Unlike the inclusion-exclusion proof of p2] , we give a simple analytic proof. 

Theorem 3 Let p(z) be any monic polynomial of degree n and let A4 n be the family of 
all n x n matrices over {—1, +1}. Then 



where PerX is the permanent of X and 7r(X) is the product of the entries of X . 

Proof. If X = (xij) is an n x n matrix of indeterminates, then by definition PerX = 
J2aes n To- where S n is the symmetric group and T a = £i CT (i) X2<r(2) • • •aW(n)- If t ne polyno- 
mial p(PerX) is expanded in terms of monomials in the x^j, then the only monomial in- 
volving every comes from products T ai T a2 ■ ■ ■ T an where the permutations a%, a-i, ■ ■ ■ > 
are the rows of a Latin square. That is, the coefficient of the only monomial with each 
Xij having odd degree is the number of Latin squares. Multiplying by 7r(X) turns the 
required monomial into the only one that has even degree in each x^. Now summing 
over X G A4 n causes this monomial to be multiplied by \A4 n \ = 2™ 2 while all the other 
monomials cancel out. □ 

6 Extremal graphs with respect to m(B) 

In our computations we learned the values of m(B) for each graph B G B(k, n) for n ^ 11. 
In Table H] we record the maximum and minimum values, and the number of graphs (in 
the column headed ) that achieve the minimum. The maximum is achieved uniquely 
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in all cases. Of course, for k ^ 1 the result is trivial and when k ^ n — 1 the unique graph 
has m(B) = i? n , so we omit these cases. 

In most cases, the graphs maximizing m(B) are the same as those with the maximum 
number of perfect matchings, as listed in [13] . The only exceptions are as follows, where 
the notation is that used in [IB"] : 

• For n = 7, k = 5 the graph maximising m(B) is 2 J2 © D3, 

• For n = 9, k = 6 the graph maximising m(B) is 3 J3; 

• For n = 10, k = 4 the graph maximising m(B) is J4 © 3J2; 

• For n = 11, k = 4 the graph maximising m(B) is J4 © J3 © D±. 

In the first of these cases the cited graph does, according to [IB], maximise the number 
of perfect matchings, but does not do so uniquely. 

7 Proportion of Latin squares with symmetry 

In this section we prove that the proportion of order n Latin squares which have a non- 
trivial symmetry tends very quickly to zero as n — > 00. 

Theorem 4 The proportion of Latin squares of order n which have a non-trivial au- 
toparatopy group is no more than 

n -3n 2 /8+o(n 2 )_ 

Proof. Suppose that a Latin square L = (£ij) of order n has a non-trivial autoparatopy 
group. Then by Lemma 4 in [IT], L has a autoparatopism a which fixes (pointwise) no 
more than one quarter of the triples of L. 

The number of possibilities for a is less than 6n! 3 = o(n 3n ). Given a, we can construct 
each possible L row by row. Each entry is determined either by a and a previous entry, 
or can be chosen in at most n ways. The latter possibility occurs once per orbit of a, and 
since a fixes at most \ of the triples of L, the number of orbits is at most {\ + \-\)n 2 = |n 2 . 
In total we find that there are most 

( n 3n )n 5n 2 /8 

Latin squares with non-trivial autoparatopy group. Our result now follows immediately 
from the well known lower bound for L n (see, for example, Thm 17.2 in [19]) that says 
that 

L n ^ (n\) 2n n- n2 ^ n n2 -°( n2 \ 

□ 
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Table 4: Minimum and Maximum values of m(B) 



9 



As a corollary to this last result, we see that the proportion of main classes or of 
isotopy classes or of isomorphism classes whose members have a non-trivial autoparatopy 
group is also bounded by ([I]). This is because each such class has somewhere between 1 
and 6n! 3 = o(n 3n ) members. 

Another corollary is that the number of isomorphism classes, isotopy classes and main 
classes of Latin squares of order n will be asymptotic to L n /n\, L n /n\ 3 and L n /(6n! 3 ), 
respectively. 
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